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ABSTRACT 

This  paper  discusses  the  dispersal  of  two  interacting  biological  species. 
The  dispersal  -  a  response  to  population  pressure  alone  -  is  modelled  by  the 
degenerate  parabolic  system 

ut  -  [u(u+v)x]x 
Vt  *  k[v(u+v)x]x 

in  conjunction  with  an  initial  prescription  of  the  individual  densities  u 
and  v  together  with  standard  zero-flux  boundary  conditions.  We  demonstrate 
here  the  following  interesting  feature  of  this  model:  segregated  Initial  data 
give  rise  to  solutions  which  are  segregated  for  all  time. 
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SIGNIFICANCE  AND  EXPLANATION 


t  /l  /  ^  ^0  '-  *• 

'-fferconsider^ber^a  mathematical  model  for  interacting  biological  species 

/7,-_  &-w  '■>  ,■>  <  /" 

that  disperse  as  a  response  to  population  pressure.  We'  demonstrate  here  an 
interesting  feature  of  the  model:  species  which  are  initially  segregated 


remain  segregated  for  all  time. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MAC,  and  not  with  the  authors  of  this  report. 


ON  INTERACTING  POPULATIONS  THAT  DISPERSE  TO  AVOID  CROWDING: 
PRESERVATION  OF  SEGREGATION 


*  *  *  ***  * 

M.  Bertsch  ,  M.  E.  Gurtin  ,  D.  Hilhorat  ,  and  L.  A.  Peletier 

1.  Introduction 

Consider  two  interacting  biological  species  with  populations 
sufficiently  dense  that  a  continuum  theory  is  applicable,  and 
assume  that  the  species  are  undergoing  dispersal  on  a  time  scale 
sufficiently  small  that  births  and  deaths  are  negligible. 

Granted  these  assumptions,  conservation  of  population  requires 
that 


ut  -  -div(uq) , 

(1.1] 

vt  »  -div(vw) , 

where  u(x,t)  and  v(x,t)  are  the  spatial  densities  of  the 
species,  while  the  vector  fields  q(x,t)  and  w(x,t)  are  the 
corresponding  dispersal  velocities. 

We  restrict  our  attention  to  situations  in  which  dispersal 
is  a  response  to  population  pressure  and  express  this  mathematically 
by  requiring  that  the  dispersal  of  each  of  the  species  be  driven 
by  the  gradient  V(u+v)  of  the  total  population,1  u+v.  We 
therefore  assume  that 
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q  =  -k^fu+v) 


(1. 


w  =  -k2V (u+v) , 


with  strictly-positive1  constants,  and  this  leads  to 

2 

the  system 


ufc  =  k^divluVIu+v) ] , 
v  =  k2div(v7(u+v) ] . 


a. 


For  convenience,  we  limit  our  attention  to  one  space-dimension 
and  we  choose  the  time-scale  so  that  =  1.  Then  writing  k  =  k2 
we  have  the  system 


ut  =  lu(u+v)x]x' 
vfc  =  k[v(u+v)xlx. 

We  shall  suppose  that  the  two  species  live  in  a  finite  habitat 


(1. 


n  =  (-L,  L)  ,  L  >  0; 


that  individuals  are  unable  to  cross  the  boundary  of  n, 

u(u+v)x  =  v(u+v)x  =  0  for  x  =  ±L,  t  >  0;  (1. 

and  that  the  two  populations  are  prescribed  initially. 


XThe  system  (1.3)  with  k2  =  0  was  studied  by  Bertsch  and 
Hilhorst  [3]  and  by  Bertsch,  Gurtin,  Hilhorst,  and  Peletier  [4]. 

2Gurtin  and  Pipkin  15] .  See  also  Busenberg  and  Travis  [6].  An 
alternative  theory  was  developed  by  Shigesada,  Kawasaki,  and 
Teramoto  [7] .  This  theory  is  discussed  in  [4]  and  [5] . 


u (x, 0)  =u0(x),  v(x, 0)  *  v0(x>  for  x  €  0. 


In  this  paper  we  shall  study  the  problem  (1.4) -(1.6) 
for  initial  data  which  are  segregated  in  the  sense  that,  for 
some  a  fc  fi  < 

Uq (x)  s  0  for  x  >  a,  vQ(x)  *  0  for  x  <  a. 

As  our  main  result  we  establish  the  existence  of  a  solution 
in  which  the  two  species  are  segregated  for  all  time.  This 
result  is  quite  surprising1  as  it  is  independent  of  the 
initial  distributions  of  the  species  and  of  the  ratio  k  of 
their  dispersivities. 


^■Actually,  Gurtin  and  Pipkin  [5]  gave  a  particular  solution  to 
(1.2)  -  corresponding  to  initial  Dirac  distributions  -  in  which 
the  two  species  are  segregated  for  all  time.  Being  a  specific 
solution,  it  is  not  clear  from  this  result  whether  "preservation 
of  segregation"  is  a  generic  property  of  the  equations  (1.4). 

2 

Granted  they  are  segregated. 


2.  The  problem.  Results. 


We  shall  use  the  notation 


IR+  =  (0,oo)  ,  q  =  n  x  ir+,  qt  =  a  x  ( 0 ,T)  , 


and,  for  any  function  f:  Q  -*  XR, 

Q+(f)  =  interior [ (x, t)  6  Q:  f(x,t)  >  0]. 


Our  problem  consists  in  finding  functions  u(x,t)  and  v(x,t) 
on  Q  such  that 


ut  =  tu(u+v)x3x 


vfc  =  k[v(u+v)x3x 


in  Q, 


( 2 .  i; 


(i) 


u (u+v) x  =  v (u+v) x  =0  on  dO  X  JR+, 


^u(x,0)  =  Uq (x)  ,  v (x , 0 )  =  v q (x)  in  0. 


(2.2; 

(2.3! 


We  shall  assume  throughout  that: 

(Al)  k  >  0,  u0,vQ  2  u0'v0  ^  C(I3)j 

(A2)  the  initial  data  are  segregated,  so  that  (1.7)  holds 
for  some  a  6  Q; 

(A3)1  each  of  the  sets  [x:  uQ(x)  >  0]  and  [x:  Vq(x)  >  OJ 
is  connected. 

The  purpose  of  this  paper  is  to  establish  -  for  such  segregated 
initial  data  -  solutions  of  (I)  which  are  segregated  for  all  time. 

1We  make  this  assumption  for  convenience  only. 


Proceeding  formally,  let  (u,v)  be  a  segregated  solution.  Then 
the  sets  Q+(u)  and  Q+(v)  are  disjoint;  hence  (assuming 
u , v  0) 

u  *  0  in  Q+(v),  v  a  0  in  Q+(u), 

and,  by  (2.1)  , 

ufc  *  <uux>x  in  Q+(u>»  vt  ■  Mwx) x  in  Q+(v).  (2.4) 

Thus  where  positive  u  and  v  obey  porous-media  equations. 

As  is  well  known,1  solutions  of  the  porous-media  equation  may  not 
be  smooth,  and  for  that  reason  it  is  advantageous  to  work  with  a 
weak  formulation  of  Problem  (I).  This  is  reinforced  by  the 
observation  that  (I)  is  a  free-boundary  problem  and  conditions 
at  the  free  boundary  are  generally  indigenous  to  a  weak  formulation, 
not  required  as  separate  restrictions.  (The  free  boundary  is  the 
set 

3  =  ldQ+(u)  U  dQ+ (v) )  0  Q  (2.5) 

which  separates  the  region  with  u(x,t)  >  0  from  that  with 
u(x,t)  *  0  and  separates  the  region  with  v(x,t)  >0  from 
that  with  v(x,t)  =  0.) 

With  this  in  mind,  assume  for  the  moment  that  (u,v)  is  a 
smooth  solution  of  (I) .  If  we  multiply  (2.1)  by  an  arbitrary 
smooth  function  )(x,t),  integrate  over  QT,  and  use  (2.2) 

^f.,  e.g.,  the  survey  article  of  Peletier  [8]. 


I 


and  (2.3) ,  we  arrive  at  the  relations 


J*  lu(T)iMT)  -  Un*(0)  ]  «  J*  Cu*  -  u(u+v)  *  J, 
n  u  qt 


lv(T)*(T)  -  v0*(0)3  =  J  [v*t  -  kv(u+v)x*x3, 

qt 


(2.( 


where  we  have  used  the  notation  u(t)  *  u(*,t),  etc.  We  shall 
use  (2.6)  as  the  basis  of  our  definition  of  a  weak  solution. 


Definition.  A  (weak)  solution  of  Problem  (I)  is  a  pair 
(u,v)  with  the  following  properties: 

(i)  u,v  €  LCC(Qt)  for  T  >  0;  u(t)  ,v(t)  6  L°°(n)  for  t  ^  0; 
(U+v)2  6  L2(0,T;  H 1 ( 0)  )  for  T  >  0; 

(ii)  u(t)  ,v(t)  >  0  almost  everywhere  in  £)  for  t  >  Of 

(iii)  u  and  v  satisfy  (2.6)  for  all  *  €  C1 (Q)  and  T  >  0 

If,  in  addition,  there  is  a  continuous  function  [0,c»)  -♦  n 

such  that,  given  any  t  >  0, 


v  (x ,  t)  =  0 

for 

-L  <  x  <  |(t)  , 

u  (x ,  t)  =  0 

for 

C(t)  <  X  <  L, 

then  (u,v)  is  segregated .  We  will  refer  to  §  as  a  separation 


curve . 


Remarks : 

1.  The  terms  u(u+v)x  and  v(u+v)x  in  (2.6)  are  defined 


as  follows: 


u(u+v)x  = 


5  35  i£  "  >  0 


if  u  *  0 


and  similarly  for  v(u+v)x.  Then,  since  u/(u+v)  <  1,  while 
t(u+v).2]x  €  L2(Qt),  we  have 


u(u+v)x  €  L  (Qt)  for  T  >  0. 


2.  The  integral  identities  (2.6)  imply  that,  as  t  -♦  0, 


u ( t)  -*  un,  v(t)  -*  vft  weakly  in  L  (ft)  ; 


i.e.,  that 


[u(t)-uQ]i|i  -•  0,  £  lv(t)-vQ)*  -»  0 


for  all  f  €  L  ift) .  To  verify  (2.8)  we  simply  apply  (2.6)  with 

4  €  C1  (?3)  (independent  of  time)  .  This  yields  (2.8)  for  t|i  6  C1(?5) 

2 

and  hence  -  using  a  standard  argument  -  for  i)  6  L  (ii)  . 

3.  The  choice  ii  s  1  in  (2.6)  leads  to  the  global 
conservation  laws 


J  u(t)  *  ;  u  ,  J  v(t)  =  J  v 

n  n  0  n  ci  0 


for  t  >  0 . 


*.v. V-/.  - 


Remark .  It  is  important  to  emphasize  that  we  have  established 
uniqueness  only  within  the  class  of  segregated  solutions.  Thus 
we  have  not  ruled  out  the  possibility  -  for  segregated  initial 
data  -  of  solutions  which  mix.  We  conjecture  that  this  cannot 
happen . 


Theorem  2.  Let  (u,v)  be  the 


ted  solution  of 


Problem  (I) .  Then: 

(i)  u  +  v  €  C(Q); 

(ii)  ufc  =  (uux)x 

classically1 

in  Q+(u); 

(iii)  v t  =  k  (wx)  x 

classically 

in  Q+(v). 

L,  .  .  GO 

That  is,  u  is  C  on 

pointwise . 

Q+(u)  and 

there  satisfies 

Our  next  theorem  is  concerned  with  the  free  boundary  3 
(cf.  (2.5)).  In  view  of  (2.4),  the  portion  of  3  along  which 
the  two  species  are  not  in  contact  should  have  properties 
similar  to  those  of  the  free  boundary  for  the  porous-media  problem 

fpt  =  (ppx>x  in  Q' 

(PM)  J  pp  =0  on  S!)  *  1+, 

X 

v_P  ( x ,  0 )  =  pQ(x)  i°  ft- 

As  is  known, ^  when  the  initial  data  have  the  form 

P q  (x)  >0  in  (a1(a2),  PQ  ( sc)  =  0  otherwise, 

-L  <  a^  <  a2  <  L,  the  free  boundary  Q  H  SQ+(p)  consists  of 
two  continuous,  time-parametrizable  curves,  one  emanating  from 
a1#  one  from  a2>  If  b(t),  0  <  t  <  T^,  designates  the  curve 
from  a^  (resp. ,  a2>  ,  then: 

( F ^ )  b (t)  =  b ( 0)  on  (0,Tfa)  for  some  6  (0,Tb); 

( F ^ )  t>(t)  is  C1  and  strictly  decreasing  (resp.,  strictly 

increasing)  on  ( ) ; 

( F 3 )  b(T~)  €  dn. 

This  discussion  should  motivate  the  following  definition  in  which 
"FB"  is  shorthand  for  "free  boundary". 

^Cf.,  e.g.,  the  review  article  by  Peletier  (8).  See  also 
Aronson  and  Peletier  (9). 


Figure  1.  b  is  a  left  free-boundary  curve  extending 
to  an.  £  is  an  internal  free -boundary 
curve  that  is  right  up  to  time  T  . 


Definition.  An  P*  curve  is  a  continuous  function 

b:  [0,tb)  -»  0 

(t^  may  be  co)  .  Moreover: 

(i)  b  is  internal  if  tb  *  ocj 

(ii)  b  is  left  (resp.,  right)  up  to  time  Tfa  €  [0,^] 
if  (FjMFj)  hold; 

(iii)  b  extends  to  an  if  b  is  right  or  left  up  to  time 
tfa  with  tb  <  co,  and  b(tb)  €  SO. 

Let  b:  (0,tfa)  -*  n  be  an  FB  curve  and  let  g:  Q  -*  JR.  Then 
FB  conditions  with  velocity  q  hold  from  the  left  (resp. ,  right 
on  b  if  given  any  t  €  (0,tb)  at  which  b  is  C1, 

b'  (t)  =  g(b(t)+,  t)  (resp.,  b’(t)  *  g(b(t)~,  t)). 

Theorem  3.  Let  (u,v)  be  the  segregated  solution  of  Problem 
(I) .  Then  there  exist  FB  curves  bu'bv»£u'£v  with  the 
following  properties: 

( i)  1  bu  <  £u  <  £y  <  bv  with  by  and  £y  forming  the 

boundary  of  Q+(u)  in  Q,  by  and  £y  forming  the  boundary 

of  Q+(v)  in  Q; 

(ii)  bu  and  bv  extend  to  an,  with  bu  left  and  by 

right; 


Here  each  inequality  is  assumed  to  hold  at  those  times  at'  which 
the  underlying  functions  are  defined. 


(iii)  Cu  and  Cy  are  internal,  and  there  is  a  time 
T  €  [0,os)  such  that  £u  is  right  up  to  time  T,  £y  is 
left  up  to  time  T,  and 

Cu<t>  *  Cy(t)  C(t)  on  [T ,co) 

with  £  €  C1^,^);  in  addition. 

(u+v )  ( C ( t )  , t)  >0  for  t  >  T;  (2.10) 

(iv)  F3  conditions  with  velocity  -ux  hold  from  the 
right  on  bu,  from  the  left  on  £u; 

(v)  FB  conditions  with  velocity  -kvx  hold  from  the 
right  on  Cy,  from  the  left  on  by. 

Remarks. 

1.  The  curve  £  marks  that  portion  of  the  free  boundary  on 
which  the  two  species  are  in  contact.  By  (2.10),  the  functions 
u  and  v  suffer  jump  discontinuities  across  £;  more  precisely 
for  t  >  T, 

u(C(t)',t)  >  0,  u(£(t)+,t)  -  0, 
v(£(t)-,t)  =  0,  v ( £ ( t) + , t)  >  0. 

Further,  (iii)-(v)  of  Theorem  3  in  conjunction  with  the  continuity 
of  u+v  imply  that,  for  t  >  T, 

u  (  £  ( t)  ” , 
ux(£(t)_,  t) 

-12- 


t)  =  v  (  £  ( t)  +  ,  t)  , 

=  kvx(£(t)+,  t)  =  -£' (t) .  (2.11) 


2.  The  results  (iv)  and  (v)  assert  that  each  of  the  "fronts" 
bu (t) »  Cu(t),  Cy(t),  by(t)  propagates  with  the  velocity  of 
individuals  situated  on  it;  condition  (2.11)  is  the  requirement 
that  at  the  contact  front  C(t)  the  two  species  move  together. 

Our  final  result  concerns  the  asymptotic  behavior  of  segregated 
solutions.  Proceeding  formally,  suppose  that  ^(x)  ,  v^x)  is  an 
equilibrium  solution  of  Problem  (I).  Then  (2.1)  and  the  boundary 
conditions  (2.2)  yield 

H»(Hx>+vco>’  ’  voo('to+voo)’  =  0  in 

hence 

( (\3tx.+  v^  )  2y  *  0  in  n 

and 

+  vaj  *  constant. 

If  Hj,  and  v^  are  segregated  with  habitats  in  [-L, 
and  [x^.,1,],  respectively,  then  there  exists  a  constant  p 


such  that 


\ 

rp 

if 

x  6  (-L,  j^j) 

Wx»  =  < 

ii  ? 

if 

x  *  ^0 

Lo 

if 

x  €  (^J.,  L) 

1 

r° 

if 

x  €  (-L,  ^ 

voo‘x)  -  H 

if 

x  -  ^ 

1 

L  ? 

if 

x  €  (XL- ,  L)  . 

Moreover,  if  the  equilibrium 

solution  (tfec.v^) 

is  reached  from 

the  initial  data  (u.,v0). 

the 

conservation  law 

(2.9)  implies 

that 


P 


U  +  V 
2L  ' 


*=c 


p 


where 

That  these  formal  calculations  are  indeed  correct  is  a 
consequence  of  the  following  theorem. 


(2.13) 


(2.14) 


Theorem  4.  Let  (u,v)  be  the  segregated  solution  of  (I) . 
Then,  as  t  -♦  co , 

i(t)  -*  u(t)  ■*  v(t)  -*  v^., 

the  latter  two  limits  being  pointwise  in  nNi^.3.  Here 
C(t)  is  the  contact  front  (cf.  Theorem  3),  while  Xj.,  , 
and  v^  are  defined  in  (2 .12 )  ,  ( 2 .13 ) . 


-16- 


Reformulation  of  the  problem. 


Let  (u,v)  be  a  sufficiently  smooth  segregated  solution  of 
Problem  (I)  and  define 


2(x,t)  -  -U  +  f  (u(y,t)  +  v(y,t) ]dy; 
-L 


z(x,t)  represents  the  total  population,  at  time  t,  in  the 
interval  [-L,x] .  In  view  of  the  conservation  laws  (2.9), 

z(-L,t)  -  -U,  2 ( 5 (t) ,  t)  -0,  z(L,t)  -V,  (3.2) 

where  U  and  V  are  defined  in  (2.14),  so  that  separation  curves 
§(t)  (cf.  (2.7))  are  level  curves  z(S(t),  t)  -  0;  in  fact, 

2 (x , t)  <  0  for  x  <  5(t) , 

<3.3  ] 

z(x,t)  >0  for  x  >  5 (t) . 

Further,  if  we  differentiate  (3.1)  with  respect  to  t  and  use 
(2.1),  (2.2),  and  (2.7),  we  find  that 


Thus  defining  c:  IR  IR  by 


s  <  0 


s  >  0 , 


we  may  use  (3.3)  to  reduce  (3.4)  to  the  single  equation 


c(z)t  "  Vxx 


on  all  of  Q.  Therefore,  if  we  write 


z0(x)  *  -U  +  J‘  (uQ+v0), 
■“L 


we  are  led  to  the  following  problem  for  z(x,t): 


fc(z)t  “  zxzxx  in  Q' 

(IQ  <z(-L,t)  =  -4  z  (L,  t)  -V,  t  >  0, 


^z(x,0)  *  z0 (x)  in  n. 


We  assume,  for  the  remainder  of  the  section,  that  c 

and  zQ  are  defined  by  (3.5)  and  (3.6)  ,  and  that  (A1)-(A3) 

are  satisfied. 


Problem  fH)  under  hypotheses  more  general  than  ours,  has 
been  analyzed  in  [10].  We  shall  simply  state,  without  proof, 
a  version  of  the  results  of  [10]  appropriate  for  our  use.  with 
this  in  mind,  we  first  define  what  we  mean  by  a  solution;  in  that 
definition,  and  in  what  follows,  ^(x)  designates  the  unique 
equilibrium  solution  of  (H)  : 


y*1 


(U+V) (x+L)  _  „ 
2L 
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Definition.  A  (weak)  solution  of  Problem  (I)  is  a  function 
z  €  C([0,oo);  w  1,0°(fl))  with  the  following  properties: 

(i)  z(*,t)  -  z^CHjcn)  for  t  >  0; 

(ii)  zt  6  L2(Qt)  for  T  >  0; 

(iii)  for  all  *  6  C1  (Q)  with  *  -  0  on  &n  x  (0,oo) 
and  all  T  >  0, 

S  Ic(z(T))*(T)  -  c(zo)*(0)3  -  J*  £c  (z)  *.  -  i(z  )2*  J. 

0  Qt  t  *  x 


Theorem  5  ( (10] )  .  Problem  (U  has  exactly  one  solution  z. 
Moreover : 

(i)  zx  €  C(Q)  with  zx  >  0/ 

(ii)  c(z)t  ■  zxzvx  classically  in  Q+(zx)? 

(iii)  Q+(zx)  is  the  union  of  the  sets 

Q2  :=  [(x,t)  €  Q:  -V  <  z(x,t)  <  O], 

Q2  {(x,t)  6  Q:  0  <  z(x,t)  <  Vj, 

and  there  exist  free-boundary  curves  b^.bjrC^r^  such  that 
(i)  -  (v)  of  Theorem  3  hold  with  Q+(u) ,Q+(v) ,bu,bv,Cu,Cy 
replaced  b^  Qi'^2,bl'b2' resPectivelv ■  with  u+v  in  (iii) 
replaced  by  z  ,  and  with  u  and  v„  in  (iv)  and  (v)  replaced 


(v)  given  any  xQ  €  n  with  z^(xQ)  >  0,  there  exists 
a  continuous  function  §:  (0,oo)  -*  Q  such  that 

t(x,t)  €  Q:  z(x,t)  *  z0(xq)3  -  Ux,t)  €  Q:  x  -  |(t)  )} 

moreover ,  5  €  c1(0,oo)  and,  for  t  >  0, 

5' (t)  -  -K2xx(5(t) ,  t) , 

where  K  =  1  or  k  according  as  Zq(xq)  <  0  or  zQ(x0)  >  0. 

The  next  result  asserts  the  equivalence  of  Problems  1  and 
(I)  and,  when  combined  with  Theorem  5, 

yields  the  validity  of  Theorems  1-4. 

Theorem  6.  Problems  (I)  and  (jl)  are  equivalent; 

(i)  Let  z  be  a  solution  of  Problem  (2)  and  define  u 
and  v  on  Q  by 

u (x , t)  =  zx(x,t) ,  v { x , t)  =  0  if  z(x,t)  <  0, 
u (x , t)  =  0,  v ( x , t)  *  zx(x,t)  if  z(x,t)  >  0,  (3.8) 

u (x , t)  =  v ( x , t)  =  ^  zx(x,t)  if  z (x.t)  *  0. 

Then  (u.v)  is  a  segregated  solution  of  Problem  (I) . 

(ii)  Conversely,  let  (u.v)  be  a  segregated  solution  of  Problem 


Proof. 


(i)  Let  z  be  a  solution  of  (Z)  and  define  (u,v)  through 
(3.8).  By  Theorem  5(i) ,  the  only  nontrivial  step  in  showing  that 
(u,v)  solves  (1)  is  proving  that  u  and  v  satisfy  the  integral 
identities  (2.6).  We  shall  only  verify  the  first  of  (2.6); 
the  verification  of  the  second  is  completely  analogous. 

For  convenience,  we  write  b  ■  b^,  C  *  ^  for  the  FB 
curves  established  in  Theorem  5,  and  we  extend  b(t)  continuously 
to  (0,oo)  by  defining  b(t)  *  -L  for  t  £  t^.  By  Theorem  5(iii) 
and  (3.8), 

supp  u(t)  «  (b(t) ,  C(t) ) . 

Choose  e  >  0  sufficiently  small  and  let  bE(t)  and  (t) , 
respectively,  be  the  level  curves  z  ■»  -0  +  e  and  z  «  -c 
(cf.  Theorem  5(v)).  Then,  by  Theorem  5 ( ii)  and  (3.8), 
ut  “  <uVx  classically  and  v  *  0  must  both  be  satisfied 
in  a  neighborhood  of  any  (x,t)  such  that  b£(t)  <  x  <  C£(t) 
and  t  >  0.  Further,  Theorem  5(v)  yields 

b'(t)  =  -u  (b  (t),  t),  C'(t>  =  -u  (w  (t),  t) 

c.  A  t  b  A  b 

for  t  >  0.  Thus,  choosing  6  >  0,  the  identity 

t 

J  f '  (T)dT  =  f  (t)  -  f  (6) 

6 


applied  to 


r\ 


C£(t) 


f(T)  =  J*  U(T)*(T) 


b£(T) 


yields,  when  i }  €  C  (Q) , 


C£(t) 


m  C£U)  t  p£(T|  -j 

J  u(t)i|i(t)  -  J  utiJiMS)  *  J  |  /  (u*  -u(u+v)  *  )dx>dT. 

(t)  be(6)  6lb£(T)  J  (3.9) 


£„(t) 


Next,  since  b£(t)  i  b(t)  and  C£(t)  '  C(t)  as  e  *  0  for  each 
t  €  (Q.,«2)  ,  it  follows  from  Lebesgue's  dominated  convergence  theorem 
that  (3.9)  holds  with  b£  and  C£  replaced  by  b  and  £.  Also, 
since  z x  €  C(Q)  ,  it  follows  that  zx(6)  "*  zq  in  C(I5)  as 


6  i  0  and 


C16) 

J  u  ( 6)  i)  ( 6)  -* 
b  ( 6) 


£  CO) 

V(°> 

b  0)  0 


as  6 


Thus  a  second  application  of  Lebesgue's  theorem  yields 


c<t)  c 

f  u(t)  ijl  (t)  -  , 

b  (t)  b 


H0)  tty*’  1 

J  uQi)(0)  =  J  ]  J  (u*  -  u(u+v)  *  IdxldT, 

b  ( 0)  u  0  |_b(T)  T  x  x  j 


and,  since  u(t)  =0  on  Q\(b(t) ,  £(t)),  the  first  of  (2.6) 
follows . 

(ii)  Let  (u,v)  solve  (I) .  We  are  to  prove  that  z  -  defined 

by  (3.1)  -  solves  (H)  .  Choose  T  >  0.  Then  u,v  6  L  (Q^)  and 

hence,  by  (2.9)  and  the  definition  of  ,  z(-,t)  -  z^  €  Hg(Q). 

2  2  1 

Note  also  that,  since  zx  =  u  +  v  and  (u+v)  6  L  (0 ,T;H  (fi) )  , 

it  follows  that 


•  i 
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€  l/(0,T;ir(n))  , 


(3.10) 


2  2 

where  z *  (zx>  . 


Our  next  step  will  be  to  establish  the  integral  identity  (3.7). 
Thus  choose  x  6  C* (Q)  with  X  =  0  on  50  x  (0,co)  and  take 


*  (x,t)  =■  j  x(y » t)  dy 
-L 


in  (2.6);  in  view  of  (2.7),  (3.1),  and  (3.6),  the  result  is 

lip  5(0)  t  5(t)  ,  , 

X  ZX(T)*(T)  -  /  z£*<0)  -J  S  £zx*t  '  I(2x)xxJdxdt' 

— L  — L  0  -L 


L  L  T  L  .  ~ 

J  z  (T)  ♦  (T)  -  X  Z'*(0)  -X  X  U„*t  -  fuXxJdxdt. 

§ (T)  x  §(0)  0  0  C(t)  x  z  1  x  x 


(3.11) 


Since  *x  *  x,  t(-L,t)  ■  0,  and  X(±L,t)  •  0,  while  z(x,t)  satisfies 
(3.2)  and  (3.3),  if  we  integrate  (3.11)  by  parts  and  then  add  the 
first  of  the  resulting  relations  to  k_1  times  the  second,  we  arrive 
at  (3.7)  (with  i)  replaced  by  x)  • 

2 

We  have  only  to  show  that  zfc  €  L  (QT> .  But  this  follows  from 
(3.10)  and  the  fact  that,  by  (3.7),  c(z)t  *  ^(zx>x  in  the  sense 
of  distributions  on  Q_.  This  completes  the  proof  of  Theorem  6. 


Remarks.  Open  problems.  Conjectures 


1.  Problem  (I)  with  nonsegregated  initial  data  is  open.  Here 
the  problem  does  not  reduce  to  a  free  boundary  problem  for  a  single 
scalar  field  z,  as  one  must  solve  the  system  (2.1)  in  regions  of 
interaction  (cf .  Remark  2) . 

2.  The  system  (2.1)  with  k  =  1  is  far  simpler  to  analyze. 
There  the  total  density  p  =  u  +  v  satisfies  (PM)  with  initial 
data  _P0  =  uQ  +  vQ,  and  once  p  is  known  (2.1)  are  linear 
hyperbolic  equations  for  u  and  v: 


ut  =  (upx)x' 


vt  =  (vpx> 


X  X 


(cf.  [5]).  Using  this  reduction  one  can  prove  uniqueness  within 
the  class  of  all  solutions  (as  opposed  to  all  segregated  solutions),  and 
one  can  show  that  solutions  which  begin  mixed  remain  mixed  for  all 
time,  including  t  =oo.  (Details  will  appear  elsewhere.) 

3.  Assume,  in  place  of  (1.2),  that  the  dispersal  of  each 
of  the  species  is  driven  by  a  weighted  sum  of  the  densities;  i.e., 
that  (in  one  space-dimension) , 


q  =  -<kuU  +  k12v)x 


w  =  -<k21u  +  k22v)x 


(4.1) 


with  all 


.2 


This  constitutive  assumption,  when  combined  with  the  conservation 
law  (1.1)  and  corresponding  zero-flux  boundary  conditions,  leads 
to  the  problem 


(nr) 


ut  =  lu(kllu  +  k12v)x]x 

in  Q, 

lvt  *  tv(k21u  +  k22v,x]x 

u(kllU  +  k15v)v  =  v(knu  +  k  ,,v)  v  =  0  on  dC)  x  m+, 


(4.3) 


22  '  x 


u(x,0)  —  Uq (x) ,  v(x,0)  ■  vQ(x)  in  n. 


This  formulation  is  greatly  simplified  if  we  define  new  independent 
variables 


a(x,t)  -k^ufx.t),  P  (x ,  t)  =  k12v(x,t) 


and  new  constants 


k  = 


22 

C12 


H  = 


k  k 

*12*21 
]r  V 

*11*22 


for  then  (m)  reduces  to 


(12) 


at  =  la(a+g)x}x 

Pt  =  k[p(pa  +  0)X3X 


in  Q, 


|  +  P)x  =  P(^a  +  P)x  =  0  on  dn  x  ir  , 


(4.4) 


^(XjO)  =  <Xg(x)  ,  P(x,0)  =  Pq(x)  in  Q, 


with 


4.  Consider  Problem  (HI),  or  equivalently  (nz)  .  The  terms  on 
the  right  side  of  (4.3)  involving  second  derivatives  are 


Writing  A  =  A{u,v)  for  the  coefficient  matrix  in  (4.5)  and 
confining  our  attention  to  u  >  0,  v  >  0,  we  conclude  from  (4.2) 
that  there  are  exactly  three  possibilities  for  the  eigenvalues 
X^  <  X2  °f  A,  namely: 

(i)  >  0,  X2  >  0;  (ii)  =  0,  X2  >  0;  (iii)  X1  <  0,  X2  >  0. 

Moreover,  writing  K  for  the  matrix 

K  -  <k..>, 

it  is  not  difficult  to  verify  that 

(i)  Xj^  >  0,  X2  >  0  <=*>  det  K  >  0, 

(ii)  Xx  =  0,  X2  >  0  <=$>  det  K  =  0, 

(iii)  X^  <  0 ,  X2  >  0  <- — >  det  K  <  0. 

We  consider  the  three  cases  separately. 

Case  (i)  (det  K  >  0) .  Here  the  system  (4.3)  is  degenerate  parabolic 
as  it  is  parabolic  when  u  >  0  and  v  >  0,  but  not  when  uv  =  0. 
Because  of  this  property,  we  expect  that  initially-segregated 
solutions  will  eventually  mix.  We  also  expect  them  to  mix  for 
another  reason.  Indeed,  assume  to  the  contrary  that  Problem  (HI) 


has  a  segregated  solution  (u,v) .  For  such  a  solution  we  would 
expect  the  two  populations  to  spread  until  they  meet,  and  then 


to  remain  in  contact  along  a  contact  front  £(t)  (cf.  Theorem  3 
and  Remark  2  following  it).  From  (4.2)  one  might  expect  that 
both  k^u  +  ki2v  an<*  ^21u  +  *22v  wou^  *>e  continuous  across 
£,  and  hence  both  zero  along  £,  a  condition  which  cannot  generally 
be  satisfied  (cf.  Remark  1  following  Theorem  3).  We  are  therefore 
led  to  the  following  conjectures  for  det  K  >  0  there  are  no 
segregated  solutions  of  Problem  (HI)  .  In  this  regard  it  would  be 
interesting  to  look  at  (HI)  with1 


e  >  0;  in  particular,  the  limit  e  ■*  0. 

Finally,  within  the  context  of  the  biological  model,  the 
off-diagonal  elements  of  K  drive  the  segregation  of  the  species, 
while  the  diagonal  elements,  by  themselves,  result  in  the  usual 
diffusive  behavior.  Since  det  K  >  0  yields  k^k22  >  *12*21' 
it  would  seem  reasonable  that  in  this  case  the  two  species 
ultimately  mix. 

Case  (ii)  (det  K  =  0)  .  Here  *  1  and  Problem  (E7)  is  identical 

to  Problem  (I) .  Thus  all  of  our  results  generalize  trivially  to 
populations  whose  interaction  is  described  by  (4.1)  with  K 
singular. 

1This  choice  of  K  arose  in  discussions  with  R.  Rostamian. 


For  the  case  det  K  *  0  we  would  like  to  call  the  system  (4.3) 


degenerate  parabolic-hyperbolic .  Indeed,  if  we  set  w  *  u  +  &, 
then,  assuming  k  >  1,  (4.4)  can  be  written  as 


wfc  -  Uw  +  (k-1)  &]wxlx, 


k(wxP)x 


i.e.,  as  a  system  composed  of  a  degenerate-parabolic  equation  and 
a  hyperbolic  equation.  The  presence  of  this  last  equation  makes 
the  discontinuity  of  u  and  v  at  the  contact  front  less 
surprising. 

In  this  case  one  can  speak  of  "passive  segregation” :  if 
the  species  start  segregated,  they  may  remain  segregated,  as  we 
have  seen  in  the  previous  sections,  and  if  they  start  mixed, 
then,  when  k  *  1,  they  remain  mixed  for  all  t  >  0  (see 
Remark  2  of  this  section) . 


Case  (iii)  (det  K  <  0) .  The  system  (4.3)  is  now  not  parabolic, 
and  Problem  (II)  is  probably  not  well  posed.  Since  the  off-diagonal 
terms  in  K  dominate  in  this  case ,  one  might  expect  a  tendency 
towards  segregation,  even  in  a  mixed  population. 


5.  The  system  (2.1)  with  k  -  0  was  studied  in  [3]  and  [4). 
There  v(x)  ■  vQ (x)  and  the  problem  reduces  to  solving  (2.1) 

(2.2) and  (2.3) In  this  case,  even  with  segregated  initial 
data,  solutions  eventually  mix,  an  apparent  contradiction  in 
behavior.  The  limit  k  -*  0  in  Problem  (I)  would  therefore  be 
interesting. 
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ABSTRACT  (continued) 


in  conjunction  with  an  initial  prescription  of  the  individual  densities  u 
and  v  together  with  standard  zero-flux  boundary  conditions.  We  demonstrate 
here  the  following  interesting  feature  of  this  model:  segregated  initial  data 
give  rise  to  solutions  which  are  segregated  for  all  time. 
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